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IDENTIFICATION  ANO  SIGNIFICANCE  OF  RESEARCH  OBJECTIVE 


I  . 


The  research  objective  is  to  theoretically  assess 
engineering  feasibility  of  manufactur i ng  low  cost,  low  loss 
millimeter  waveguide  systems  ( 10  mm  -  1  mm)  using  methods  from 
fiber  optics  research  developed  for  production  of  glass  fiber 
lightguides.  To  more  completely  utilize  such  optical  fiber 
construction  concepts,  the  initial  research  effort  was  a 
theoretical  study  of  graded  index  of  refraction  Solid  Guides  of 
circular  cross  section  and  0-type  Guides  [1,2]  of  circular  cross 
section  with  a  graded  index  of  refraction  wall.  More  complex 
wave  guiding  structures  such  as  non-circular  cross  sections 
(e.g.,  elliptical),  could  be  manufactured  using  these  techniques. 

A  number  of  significant  technologies  can  be  implemented  upon 
the  successful  development  of  low  cost,  low  loss  millimeter 
waveguide  systems.  These  technologies  include  wide  bandwidth 
commun i cat  ions  links  which  can  be  EMP  hardened;  sma 1 1 , 1 i ghtwe i ght 
radars  for  autonomously  guided  weapons  and  area  survei I  lance 
using  imaging-mode  operation;  sensors  for  robotics  application; 
and  free  space,  line  of  sight  communications  systems  (e.g.,  man 
to  man,  aircraft  to  aircraft,  ship  to  ship,  etc.)  offering 
security  and  anti  jam  capability.  A  further  significant  advantage 


for  such  low  cost,  low  loss  millimeter  waveguide  systems  is  the 
possibility  of  developing  analog  signal  processing  devices  of 
reasonable  size,  including  periodic  structure  [3]  and  transversal 
filters  [4],  In  fact,  the  original  realization  of  the  transversal 
filter  was  implemented  with  coaxial  cable.  The  implication  is 
that  systems  operating  at  these  wavelengths  with  low  loss 
waveguides  can  achieve  signal  processing  without  conversion  to  an 
intermediate  frequency. 


2 


r 


II.  TECHNICAL  SUMMARY 

For  millimeter  wavelength  systems  ( l 0  mm  -  1.0  mm),  guides 
are  presently  fabricated  to  high  mechanical  precision,  usually 
out  of  costly  materials  such  as  coin  silver.  Costs  for  such 
guides  is  in  the  order  $500  per  meter  of  transmission  link  run  if 
the  number  of  connections  is  small;  and,  the  advantage  afforded 
by  operation  at  these  smaller  wavelengths  is  lost  by  the  large, 
poor  form  factor  assemblies  presently  used.  What  was  proposed  is 
to  adapt  the  construction  techniques  developed  by  workers  in 
fiber  optics  [5,6]  to  manufacture  low  loss  graded  index  of 
refraction  Sol  id  Guides  and  0-type  Guides  with  a  graded  index  of 
refraction  wall  suitable  for  use  at  millimeter  wavelengths.  In 
this  process  a  preform  of  the  desired  cross  section  parameter 
variation  is  manufactured  using  chemical  vapor  deposition  and  is 
drawn  to  the  desired  diameter;  for  optical  fiber  waveguides, 
dimensions  are  held  to  within  one  percent  for  fibers  five 
thousandths  of  an  inch  in  diameter.  Such  mechanical  precision 
should  also  be  obtained  in  the  manufacture  of  millimeter 
waveguides  by  a  similar  process.  The  techniques  of  drawing  the 
waveguides  as  is  done  for  optical  fibers  should  lead  to 
significant  reduction  in  production  cost  while  maintaining  the 
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mechanical  precision  required.  Further,  initial  considerations 
indicate  the  possibility  of  a  significant  reduction  in  the 
propagation  attenuation.  Available  rectangular  metal  pipe 
waveguides  such  as  WR5  (RG135/U)  have  attenuations  of  the  order 
25  dB  per  meter;  in  contrast,  results  from  this  study  show  that  a 
few  tenths  of  a  dB  per  meter  attenuation  are  readily  attainable 
for  mi  1 1 imeter  wavelength  operation  in  guides  constructed  of  low 
loss  dielectric  materials  such  as  watei — free  fused  silica 
(Thermal  American  Fused  Quartz  "spectras i 1 " )  and  that  lower 
values  may  be  achieved.  It  is  noteworthy  that  1  dB  per  meter  at 
one  millimeter  wavelength  corresponds  to  a  figure  of  merit  of 
1000  wavelengths  per  dB  of  attenuation,  a  figure  of  merit 
indicating  the  suitability  of  the  transmission  medium  for 
utilization  in  realization  of  transversal  filters;  corresponding 
values  for  other  transversal  filter  realizations  media  are  1000- 
5000  for  SAW,  500  for  MSW  and  500  for  copper  guide  at  microwave 
frequencies.  Since  attenuations  of  much  less  than  IdB  per  meter 
can  be  attained  in  the  proposed  guides,  this  propagation  medium 
approaches  and  exceeds  the  figure  of  merit  of  SAW  for  transversal 


f i 1 ters . 


II!  .  OESIGN  GUIDELINES 


In  this  section  we  develop  the  basic  design  criteria  used  in 
our  design  for  the  O-type  millimeter  waveguide.  The  bas i c 
structure  under  cons i derat i on  is  shown  in  Figure  1.  We  begin  by 
developing  a  simple  transmission  line  analogy,  and  then  move  on 
to  the  design  strategy  for  the  millimeter  waveguide. 

A.  Transmi ss ion  L i ne  Ana  1 ogy 

Here  we  develop  criteria  for  maximization  of  voltage 
reflection  from  a  simple  alternating,  multi-section  transmission 
line  circuit.  The  results  are  useful  in  formulating  analogous 
criteria  for  maximization  of  electromagnetic  wave  reflection  from 
alternating  multi-layered  dielectric  structures. 


Cons i der 

the 

a  1 ternat i ng 

multi 

-sect  ion 

transmi ss i on 

1  i  ne 

circuit  shown 

i  n 

Figure  2a. 

It 

cons i sts 

of  a  source 

with 

internal  impedance  Z s  connected  to  a  load  of  impedance  Z^  via  N 
transmission  line  segments.  The  transmission  line  segments  have 
alternating  character i st i c  impedances  ZQ1  and  ZQ2,  and  are  of 
alternating  lengths  X j /4  and  where  Xj  and  X2  are  the 
wavelengths  in  the  alternating  sections. 
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Using  elementary  transmission  line  analysis  it  can  be  shown 


that  the  input  impedance  seen  by  the  source  is  given  by 


so  that  if  Zqj  *  ^02  anc*  N  's  ,ar9e  we  have  either  an  extremely 
large  or  extremely  small  input  impedance.  Now  the  voltage 
reflection  coefficient  at  the  source  is  given  by 


Fi  - 


Zi  -  Zs 
2i  +  Zs 


which  has  a  magnitude  very  close  to  one  if  N  is  large. 

By  this  analysis  we  see  that  the  voltage  reflection 
coefficient  can  be  maximized  by:  1)  increasing  the  contrast 
between  the  character i st i c  impedances  ZQ1  and  ZQ2 .  and  2) 
increasing  the  number  N  of  transmission  line  sections. 


B .  Qes i gn  Parameters 


In  designing  our  millimeter  waveguide  we  will  adapt  these 
rules  for  maximizing  the  reflection  of  electromagnetic  waves 
inside  of  the  guide. 
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we  w| | 1  concentrate  on  stroctures  with  alternating  -all 
composition  and  thickness.  The  -an  thicknesses  w.  I  1  be  chosen 
to  be  a  guarter  wavelength  .n  thickness  for  waves  of  freduency  f„ 

incident  on  the  -all  at  angle  of  incidence  e„.  Referring  to 


Figure  2b,  this  means  that 


t,  = 


4  6 r |  cose. 


where  AQ  =  c/F,  -  VicO’  and  bV  SneH'5  ’aW 
e.  a  sin  '(sin0Q/  erj ) 

we  will  use  as  our  standard  test  case  in  the  rest  of  this 
document  the  following  parameters: 


fft  =  1QQ  GHz,  0n  =  80  degrees 


t  ,  =  3.285, 
r  1 


e  =  1.2 c  ,  *  3.942 
er2  r  1 


N  =  16,  32,  64  layers 
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IV.  REFLECTION  FROM  LAYERED  STRUCTURES 


In  this  section  we  develop  two  methods  of  finding  the 
reflection  coefficient  of  plane  electromagnetic  waves  of 
arbitrary  frequency  and  polarization  incident  at  any  angle  on  a 
infinite  planar  layered  dielectric  structure.  We  first  restrict 
our  analysis  to  lossless  dielectrics  here;  later  the  results  will 
be  generalized  to  lossy  dielectrics.  These  results  will  be  used 
later  in  our  ana  I ys i s  of  wave  guidance  by  a  multi-layered 
dielectric  planar  waveguide.  As  we  will  show  later,  this 
behavior  of  this  waveguide  closely  approximates  the  behavior  of 
the  multi-layered  dielectric  O-type  waveguide. 

A.  Ref  1 ect i on  from  Loss  1  ess  P 1 anar  Layers 

Consider  the  situation  depicted  in  Figure  3.  An  arbitrarily 
polarized  plane  wave  of  magnitude  A  in  free-space  is  incident  on 
an  infinite  planar  mu  1 1 i - 1 ayered  dielectric  structure.  The 
structure  consists  of  N  homogeneous  lossless  dielectric  slabs  of 
permittivity  e.  and  thickness  t.,  for  i=l,...,N.  We  seek  the 
magnitude  of  the  resultant  reflected  wave  in  free-space,  denoted 
B  in  Fi gure  3 . 
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Following  Wait  [7j  we  can  define  wave  impedances  in  each 
region  as  the  ratio  of  the  tangential  electric  field  to  the 
tangential  magnetic  field.  Electromagnetic  boundary  conditions 
are  then  enforced  by  equating  wave  impedances  on  each  side  of 
each  dielectric  boundary.  Beginning  this  process  at  the  boundary 
at  y=0  and  continuing  to  the  boundary  at  y=tN  yields  the 
following  iterative  formula  for  the  overall  reflection 
coefficient,  defined  as  the  ratio  of  B  to  A: 


ZyN  ~  z0 
ZyN  +  Z0 


ZyN  "  ZN 


ZyN-lCOs(kyNtN)  +  jZNS  1  n  (  Wn  > 
ZNcos(kyNtN)  +  jZyN-lsin(kyNtN) 


Zyi  ~  Zi 


Zyj-1cos(ky.t.)  +  jZtsin(kyttt) 
Z  j  cos  ( ky 1 1 )  t  jZyj_lSin(l<yit}) 


Z0CO5(kyltl)  JZlsfn^kvlt  1  ? 

1  ZiCOS(kyltl)  +  Jz0s i n ( ky j  t , ) 
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where 
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r 


where  [T]  denotes  the  matrix 


}j(<$1+<S2)_R2ej(<S1  <S2)  2jR  ji,  sin  6 

1  I  i. 


1  -  RZ  1_  -2  jR  j  e~^  1  sin  &2  e  j  ( 3  1 +6Z  }  -RZe  ^  ( 1  -<S2  } 


Here 


R0  = 


Z1  zo 


Z1  +  zo 


R  i  - 


Z1  -  Z2 
Z1  +  Z2 


RN  - 


Z0  -  Z2 
Z0  +  Z2 


with  Z.  defined  as  above,  and  <J  {  =  k  jtj. 

In  order  to  use  this  formulation  to  find  the  overall 
reflection  coefficient  we  must  have  an  efficient  way  to  raise  the 
matrix  [T]  to  large  powers. 

Sylvester's  theorem  as  outlined  in  Collins  [8]  allows  us  to 
express  [T]  raised  to  the  N  power  in  terms  of  the  sum  of  two 
matrices  multiplied  by  powers  of  eigenvalues  of  the  original 
matrix.  This  means  that  calculation  of  [T]  to  arbitrarily  high 
powers  can  be  accomplished  by  addition  of  two  matrices,  resulting 
in  a  vast  savings  of  time  when  N  is  large. 

This  technique  of  computing  the  reflection  coefficient 
versus  angle  of  incidence  and  frequency  was  implemented  in  the 
computer  program  WMAT.  Output  of  the  program  is 


indistinguishable  from  that  provided  in  Figures  4  through  7. 


8-  Ref 1 ect ion  from  Loss  1  ess  Cv 1 i ndr I ca I  Lavers 

In  this  section  we  use  the  wave  impedance  technique 
introduced  in  the  previous  section  to  solve  for  the  reflection 
coefficient  for  cylindrical  waves  of  arbitrary  frequency  and 
polarization  incident  at  any  angle  on  a  multi-layered  cylindrical 
structure.  Once  again  we  derive  results  for  the  lossless  case, 
which  will  be  generalized  to  the  lossy  case  later  on. 

Consider  the  situation  depicted  in  Figure  8.  An  arbitrarily 
polarized  cylindrical  wave  of  amplitude  A  is  incident  at  angle  e 
on  an  i nf i n i te I y-1 ong  cylindrical  structure.  The  structure  is 
composed  of  N  concentric  homogeneous  dielectric  layers  of 
permittivity  e j  and  inner  radius  p^,  where  i=l,...,N.  We  seek 
the  magnitude  B  of  the  resultant  reflected  cylindrical  wave. 

Once  again  we  define  the  wave  impedances  in  each  layer  as 
the  ratio  of  total  tangential  electric  field  to  total  tangential 
magnetic  field,  and  enforce  boundary  conditions  by  equating  wave 
impedances  on  both  sides  of  each  boundary.  This  process  yields 
an  interative  formula  for  the  overall  reflection  coefficient.  In 
the  case  of  TM  polarization  of  the  incident  wave,  we  obtain 


H 


(2) 


n  (kpN<V  +  KN-IHn  (kpN^N> 
PN  ■  ZN  m(2). 


H’ 


n 


(,Wn>  *  KN-lHn' ’ '<kpNPN> 


K . 


h  +  lHn' V  ~  z.lHn2>'(k 


Pilpli 


l»l>  -  zl»lHn‘><koi1. 


> 


Zpi  » 


H 


Hn2’  <  kp  i»|>  * 

(2)  , 
n 


(kpjPf)  + 


K0  » 


ZlHn2>(kplp0)  "  Z-«HlZ)'(k~.<>n> 


~pQ  n 


PlHQ' 


ZpOHil>#(kplpO>  "  2lHn,>(kp.p>0) 


z  =  2  Hn2)(kPOpO} 

0  H<2)'(kpQp0) 

where  and  H^2^  are  Hankel  functions  of  the  first  and  second 

n  n 

kind  of  order  n,  and  primes  denote  derivatives  with  respect  to 


argument . 

The 

Z.  = 

i 

k 

Pi 

u  c . 

i 

and  the 

order 

k  .  =  kn 
p  i  0 


V  Er i 


s i n20f 


and  the  order  n  denotes  the  angular  variation  e 


jn* 


of  the 


incident  and  reflected  fields  (usually  taken  as  zero). 
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In  the  case  of  a  TE  polarized  Incident  field 


a  similar 


derivation  yields 


r  = 


Z0Hn2)'(kpoV  ' 


~PN  n 


oOHN ; 


2pNHn1,|k»0»N»  -  Z0Hn‘,'<koOl>N» 


ZpN  "  ZN 


Hn2>'(k.N»N’  *  KN-lHn'>'tkpNV 


Hn2><kPN*N>  +  KN-lHn1,<kpNON> 


K  i  - 


Zi  +  lHn2>'<k„i>l»,l  '  Z»lH",<k»ll«l> 


ZpiHo')(k»i*l‘>l>  *  2l*lHn 


Zp<  *  2, 


,(Z) 
n 


<  1  )  . 

■> 

(I  ) 


(•<-.  .  ,Pi  ) 


H-‘"tkpl»i)  ♦  Kt-iHA‘ "<k0i»i > 


p  i  + 1 

Cp  i ( 


Hi2)<k»l»l>  *  kl-lHn1,|kpi».1 


K0  * 


Z’H-2)'<k»l‘>0>  -  Z»0Hn2’lkPl»0> 


‘  1  n 


~pQ  n 


:oOHn1,<kpl»0>  '  ZlHn1)-<kpl»0> 

,(2>  , 


h: 


Zp0  ~  zo 


n 


(kpQp0) 


Hn2><kpO»0> 


with  the  TE  wave  impedance  defined  by 


uu. 


Z.  = 

i 


p  > 


where  k 


is  given  as  above. 


"\ 
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This  technique  of  finding  the  reflection  coefficient  versus 
angle  of  incidence  and  frequency  was  implemented  in  the  computer 
program  CYL .  Graphical  results  of  this  program  are  provided  in 
Figures  9  and  10. 

Note  that  the  curves  for  this  cylindrical  case  are  very 
close  to  the  curves  in  the  planar  case  when  the  curvature  of  the 
cylinder  is  on  the  order  of  a  wavelength  or  more.  This  result 
implies  that  we  can  use  the  simpler  planar  case  as  a  close 
approximation  to  the  actual  cylindrical  guide. 

C.  Ref 1 ect ion  from  Lossy  Lavers 

In  this  section  we  generalize  the  formulations  of  the 
previous  sections  to  obtain  reflection  coefficients  for  lossy 
dielectric  layers.  We  first  consider  the  planar  case,  and  then 
the  cylindrical  case. 

The  generalization  to  lossy  dielectrics  in  the  planar  case 
is  a  trivial  one.  All  of  the  previous  results  hold  in  the  lossy 
case  if  we  replace  the  permittivity  e  everywhere  by  the  complex 
permi tt I v i ty , 

e  =  E'  -  je" 
c  J 


1  5 


where  e'  is  the  real  permittivity  and  e"  is  the  loss  tangent  of 
the  dielectric.  The  loss  tangent  can  also  be  written  in  terms  of 
the  conductivity  of  the  dielectric  as  em  =  o/«. 

This  generalization  was  incorporated  into  the  wave  impedance 
formulation  for  the  reflection  coefficient  from  planar  layers. 
The  modified  program  LREFL  computes  the  reflection  coefficient 
versus  angle  of  incidence  and  frequency  in  the  lossy  case. 
Sample  output  from  the  program  is  provided  in  Figure  11. 

The  generalization  to  lossy  dielectrics  in  the  cylindrical 
case  i s  obta i ned  i n  exact  1 y  the  same  way  as  above .  One  added 
complication,  however,  is  the  requirement  for  computation  of 
Hankel  functions  of  complex  arguments. 

The  program  LCYL  computes  the  reflection  coefficient  versus 
angle  of  incidence  and  frequency  in  the  lossy  case.  Sample 
output  of  the  program  is  provided  in  Figure  12. 
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V.  WAVEGUIDE  ATTENUATION  CALCULATION 


In  this  section  we  use  the  results  of  the  previous  sections 
to  calculate  the  attenuation  of  electromagnetic  waves  guided  by 
layered  structures.  We  restrict  our  attention  to  the  planar 
case;  extension  to  the  cylindrical  case  is  discussed  in  the  next 
section. 

A.  Approx i mate  Loss  Ca i cu 1  at i on 

An  simple,  approximate  waveguide  loss  calculation  can  be 
obtained  as  follows. 

Consider  a  plane  wave  of  frequency  f  bouncing  at  an  angle  © 
down  the  interior  of  the  planar  layered  structure  as  shown  in 
Figure  13.  Successive  reflections  occur  at  a  distance  d  apart, 
and  at  each  reflection  the  magnitude  of  the  wave  is  decreased  by 
the  multiplicative  factor  JRJ,  which  is  the  magnitude  of  the 
reflection  coefficient.  If  we  assume  that  the  wave  is  decreasing 
with  distance  in  the  form  e  QZ,  then  we  can  solve  for  a  as 

a  =  - 1 n (  • R  J ) /d  Np/m 

Using  the  fact  that  d  =  t  tan©  and  I  Np/m  =  8.686  dB/m,  we  obtain 


1  7 


’T 


1  n  (  !  R  |  ) 

a  =  -8.686  - 

t  tane 

Sample  values  of  a  for  the  test  cases  are  provided  in  Table 

1  . 

It  should  be  noted  carefully  that  this  calculation  is  a  very 
rough  one  due  to  numerous  approx i mat i ons  being  used.  They  are: 

1.  We  have  assumed  that  the  wave  being  guided  is  a  single 
plane  wave.  A  better  approx imat f on  is  obtfaned  by  assuming  that 
the  field  in  the  guide  is  the  sum  of  several  plane  waves 
superimposed.  However,  even  this  yields  approximate  results, 
since  we  have  no  guarantee  that  the  waves  guided  by  this 
structure  are  plane  waves. 

2.  We  have  assumed  that  we  know  a-priori  the  angle  and 
frequency  of  propagation  in  the  waveguide.  It  is  entirely 
possible  that  the  dominant  mode  in  the  guide  might  have  another 
frequency  and  might  bounce  down  the  guide  at  another  angle. 

B.  Exact  0 i spers i on  Cal cu I  at i on 

A  more  exact  analysis  of  the  dispersion  of  the  planar  guide 
is  provided  here.  We  will  make  use  of  an  equivalent  guide,  and 
solve  exactly  for  the  dispersion  relation  of  this  equivalent 
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Given  the  complex  reflection  coefficient  R  for  plane  waves 
incident  on  a  homogeneous  material  at  the  angle  6,  we  can  compute 


the  permittivity  and  conductivity  of  the  material.  A  simple 
calculation  shows 


-  D  \ 


?  2 

=  sin  9  +  cos  0  Re 


1  +  R 


/I 


a  =  -mEqCOS^S  Im 


1  +  R 

where  Re  and  Im  denote  the  real  and  imaginary  parts, 
respect i ve 1 y . 

Now  consider  the  propagation  of  electromagnetic  waves  in  the 
situation  depicted  in  Figure  14.  It  can  be  shown  that  the 
propagation  constant  for  this  system  satisfies  the  following 
characteristic  equation: 


(ky0  +  ky)sin(kyOt)  =  ZjkyOkycos(kyOt) 

where 

V  ■  /k 0  -  ‘2  •  ky  *  /  cek0  -  ‘2 

with  c  the  complex  propagation  constant  for  the  waveguide. 


This  formulation  was  used  to  calculate  c,  for  the  planar 


layered  guide,  and  the  resultant  dispersion  curves  produced 
the  program  DISP  are  given  in  Figures  15  and  16. 
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VI.  REFLECTION  COEFFICIENT  BY  RAY  TRACING 


A  computer  program  based  on  a  ray  tracing  -  transmission 
line  technique  has  been  developed  to  calculate  the  reflection 
coefficient  of  a  multiple  layer,  planar,  dielectric  waveguiding 
structure.  The  computer  program  calculates  the  amplitude  and 
phase  shift  of  the  reflection  coefficient  at  the  input  wall  for 
any  number  of  dielectric  layers  as  a  function  of  incident  angle. 
The  data  generated  by  this  program  correlates  very  well  with 
waveguide  algorithms  developed  by  other  members  of  the  study 
group.  Inputs  to  the  program  include: 


*  Number  of  Layers 

*  Illuminating  Wave  Frequency 

*  Illuminating  Wave  Angle  Incidence 

*  Quarter  Wavelength  Matching  Information  For  Layer  Thickness 

Calculation 

1.  Matched  Frequency 

2.  Matched  Incident  Angle  In  The  Layer 

3.  Layer  Dielectric  Coefficient 

*  Propagation  Mode  Type 

1.  Transverse  Electric  (TE)  -  Perpendicular 
Po 1 ar i zat i on 

2.  Transverse  Magnetic  (TM)  -  Parallel 
Pol ar i zat i on 

*  Character i st i cs  of  the  Surrounding  Media 


The  program  begins  by  calculating  the  thicknesses  of  the 
dielectric  layers  for  quarter  wavelength  matching  at  a  fixed 
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incident  angle  and  fixed  frequency  for  each  layer.  After  this 
calculation  is  completed  the  thicknesses  are  stored  for  later 
use.  Snell's  law  is  used  to  calculate  for  each  incident  angle  of 
illumination  the  necessary  incident  angles  of  propagation  in  the 
layers.  The  main  computation  then  proceeds  using  standard 
transmission  line  techniques.  Starting  at  the  load  (non- 
illuminated)  side  of  the  structure  the  program  calculates  the 
complex  reflection  coefficient  of  each  layer  interface  taking 
into  account  the  character i st i cs  of  surrounding  media  (usually 
air)  at  the  first  (load)  interface-  The  formula  used  for  calcu¬ 
lating  the  complex  reflection  coefficient  depends  upon  the  propa¬ 
gation  characteristics  (transverse  electric  or  transverse 
magnetic)  of  the  illuminating  wave. 

After  the  reflection  coefficient  is  calculated  the  input 
impedance  at  the  input  side  of  the  layer  is  computed  taking  into 
account  the  correct  propagation  distance.  This  is  then  used  to 
calculate  the  reflection  coefficient  at  the  new  interface.  The 
process  is  repeated  until  the  illuminated  interface  is  arrived 
at.  The  complex  reflection  coefficient  as  a  function  of  incident 
angle  is  then  plotted. 


A  parallel  plate  waveguide  using  dielectric  walls  can  be 
analyzed  using  this  method.  Alpha  and  Beta  can  be  calculated  for 
high  values  of  reflection  coefficient.  Plots  of  alpha  and  beta 
at  a  function  of  frequency  can  then  be  obtained. 

A.  TE  and  TM  Modes 

The  propagation  mode  types  are  transverse  electric  (TE)  and 
transverse  magnetic  (TM)  depending  upon  whether  the  electric 
vector  is  perpendicular  or  parallel  to  the  plane  of  incidence. 
The  complex  reflection  coefficient  for  these  two  cases  are  given 
in  Section  IV  A  of  this  report. 

B.  Transmi ss ion  L i ne  I nput  I mpedance 

The  input  impedance  to  each  layer  is  calculated  from  the 
following  fami 1 iar  relation  from  transmission  line  theory, 

1  +  r  exp ( -  j8d  ) 

2  in  =  - — 

1  -  T  exp(-j8dn) 

where  r  is  the  complex  reflection  coefficient,  8  is  the  phase 
constant  of  the  wave,  and  dn  is  the  electrical  path  length  in  the 
I ayer . 
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C.  Sne  1  1  s  Law  For  Transmi  ss  i  on  Anq  1  es 

The  Snell  Relation  gives  the  transmission  angle  to  the  normal  at 
a  dielectric  media  boundary. 

D.  Quarter  Wavel enoth  Laver  Match i nq 

Quarter  wave  matching  in  each  layer  is  provided  by  making  the 
layer  widths  such  that  the  wave  experiences  a  quarter  wavelength 
propagation  path  at  its  transmission  angle  in  the  dielectric 
medium.  The  matching  geometry  is  shown  in  Figure  2b. 

E .  Para  I  1  el  Plate  Waveou i de  Us i nq  Layered  D i e I ectr i c  Walls 

A  parallel  plate  waveguide  using  layered  dielectric  walls 
can  be  analyzed  using  this  method.  Alpha  and  beta  can  be  arrived 
at  for  high  values  of  reflection  coefficient  by  the  following 
formulas.  Note  that  beta  is  a  function  of  geometry  (wall 


separation)  only. 


F.  FORTRAN  77  COMPUTER  PROGRAM 


A  FORTRAN  77  program  for  this  analysis  has  been  written  for 
both  the  U.T.A.  IBM  Mainframe  and  MS-FORTRAN  for  use  on  a  PC.  In 
both  cases  the  complex  reflection  coefficient  is  written  to  a 
data  file  and  then  imported  to  lotus  123  (Ver.  2)  for  plotting. 
This  is  possible  if  the  data  file  has  a  .PRN  extension.  In  the 
mainframe  case  HAYES  SMARTCOM  II  1200  BAUD  COMMUNICATIONS 
software  is  used  to  down  1  oak  the  data.  In  the  MS-FORTRAN  case 
the  program  generates  the  data  file  directly. 

The  first  point  of  investigation  into  the  performance  of 
the  multi-layered  dielectric  waveguide  ws  the  development  of  its 
mode  characteristics.  Looking  at  a  single  boundary  between 
dielectrics,  the  TE  mode  reflection  is  given  above,  where  the 
magnetic  permeability  of  the  material  is  taken  at  unity.  In  the 
TM  mode  a  similar  expression  is  presented  above. 

The  computer  analysis  program  developed  to  model  the 
reflecting  wall  sequentially  performed  an  impedance 
transformat  ion  from  the  outside  or  load  side  of  the  multiple 
layered  dielectric  structure  to  the  inner  boundary  of  the 
reflecting  wall,  and  at  that  point  determined  the  reflection 
characteristic  as  a  function  of  input  frequency  and  incident 
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angle. 

The  first  data  was  developed  from  a  16  layer  structure  with 
5  percent  contrast  between  layers,  dsigned  to  operate  at  100  GHz 
with  an  incident  angle  of  80  degrees.  The  reflection 
character i st i cs  of  the  TE  and  TM  modes  are  shown  in  Figure  17, 
demonstrat i ng  the  pre-eminence  of  the  TE  mode  for  low  grazing 
angles  (incident  angles  above  80  degrees). 

From  this  basis,  invstigation  into  the  performance  of  a 
parallel  plate  waveguide  fabricated  with  layered  dielectric  walls 
began.  In  the  case  of  a  parallel  plate  guide  with  an  air  or 
evacuated  core,  a  plate  seperation  t,  operated  in  mode  number  n, 
the  cut-off  frequency  of  the  structure  is  given  bv  nc/2t,  where  c 
is  the  speed  of  1 i ght  in  vacuum.  In  such  a  guide  the  angle  of 
incidence  of  the  wavefront  to  the  surface  wal Is  is  given  by 

=  cos  1 (nc/2ft) 

Using  these  parametric  expressions  to  derive  the  incident 
angle  for  a  given  frequency,  the  performance  of  a  waveguide 
structure  can  be  analyzed  for  a  wide  spectrum  of  input 
frequencies  given  a  defined  geometry  specified  by  the  layer 
permittivity  and  contrast,  number  of  layers,  design  frequency  and 


26 


incident  angle,  and  the  plate  seperatlon.  There  is  a  direct 
relationship  between  the  plate  seperatlon  and  the  design  incident 
angle  at  eh  design  frequency  that  configures  the  guide  to  a 
particular  operating  point.  The  results  of  such  a  wide  band 
analysis  are  presented  in  Figure  18. 

The  attenuation  and  transmission  character i s i cs  of  the 
structure  have  also  been  observed  using  the  expressions  derived 
above.  For  the  16  the  layer  case  the  resulting  attenuation  and 
transmission  curves  are  shown  in  Figures  1 9a  and  1 9b. 

It  may  be  observed  that  the  reflection  characteristics  of 
the  model  clearly  demonstrate  a  harmonic  character,  and  the  guide 
appears  to  improve  in  performance  at  harmonic  frequencies. 
Returning  to  the  reflecting  wall  analysis  to  check  this 
observation  it  was  found  that  at  low  grazing  angles  of  10  degrees 
or  less  the  performance  was  not  degraded  between  a  64  layer  100 
GHz  design  anad  a  3  2  layer  50  GHz  design,  when  oerating  both  at 
100  GHz.  This  development  shows  the  promise  of  decreased 
difficulty  in  manufacturing  and  correspond i ng 1 y  lowered 
production  costs.  These  results  are  shown  in  Figures  20  and  21. 

Follow  on  research  will  further  investigate  the 
cons i derat i ons  of  geometry,  harmonic  operation  and  materials;  as 


well  as  transition  Into  the  cylindrical  and  elliptical  waveguide 
structures . 
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vtl.  SUMMARY  AND  CONCLUSIONS 

In  the  research  reported  here  the  original  concept  of  using 
layered  dielectric  walls  to  achieve  low  loss  propagation  media 
for  millimeter  waveguide  systems  has  been  shown  to  be 
theoret 1 ca 1 1 y  f eas i b 1 e . 

The  research  completed  as  of  this  reporting  indicates  the 

f o 1 1 ow i ng  s 

1.  there  exist  low  loss.  approximately  10  percent  of 

center  frequency  bandwidth  TE  propagation  modes  for  the 
guides  as  proposed. 

2.  the  structures  tend  to  be  excessively  large  in 

transverse  dimensions  at  the  longer  wavelengths; 
however,  as  expected,  the  structures  appear  to  be  of  a 
very  desirable  size  at  the  short  wavelengths. 

3.  the  materials  data  is  difficult  to  obtain  and 

substantiate  because  of  the  propritary  nature 
associated  with  material  development  and  the  technical 
difficulty  in  making  suitable  measurements  at  the 

proposed  frequencies. 

It  is  the  conclusion  of  the  researchers  that  these  concepts 
should  be  further  examined  theoretically  and  experimentally,  but 
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most  importantly  experimentally.  Further,  it  is  the  strong 
recommendation  that  further  research  on  these  waveguiding 
structures  should  be  focused  on  the  spectrum  above  300  GHz.  This 
means  that  these  very  low  loss  waveguide  modes  correspond  to 
technically  useful  structures  at  these  higher  frequencies. 
Finally,  this  concept  is  presently,  to  the  researchers  knowledge, 
the  only  viable  waveguiding  structure  at  these  frequencies. 
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TABLE  I 

Example:  let  f  =  1 00  GHz,  d  =  2\  0  =  88° 

or  required  1F1 

1  dB/m  .980 

0  5  dB/m  .990 

0  1  dB/m  .998 
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